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Abstract 

In 1991, Edelsbrunner and Tan gave an 0(n 2 ) algorithm for finding 
the MinMax Length triangulation of a set of points in the plane. In 
this paper we resolve one of the open problems stated in that paper, 
by showing that finding a MaxMin Length triangulation is an NP- 
completc problem. The proof implies that (unless P=NP), there is no 
polynomial-time approximation algorithm that can approximate the 
problem within any polynomial factor. 

1 Introduction 

Triangulating a set of points is one of the basic problems of Computational 
Geometry: given a set P of n points in the plane, connect them by a maximal 
set A of non-crossing line segments. This implies that all bounded faces of 
the resulting planar arrangement are triangles, while the exterior face is the 
complement of the convex hull of P. 

Triangulations are computed and used in a large variety of contexts, 
e.g., in mesh generation, but also as a stepping stone for other tasks. While 
it is not hard to compute some triangulation, most of these tasks require 
triangulations with special properties that should be optimized. Examples 
include maximizing the minimum angle, minimizing the total edge weight 
or the longest edge length. 

In this paper we consider the task of computing a triangulation whose 
shortest edge is as long as possible. We show that this problem is NP- 
complete, resolving an open problem stated by Edelsbrunner and Tan in 
1991 [4J. The proof implies that (unless P=NP), there cannot be any 
polynomial-time approximation that gets with any polynomial factor. 
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Related Work. For a broad survey of triangulations in a variety 
of settings, see the book [2] by De Loeara, Rambau, and Santos. Maxi- 
mizing the minimum angle in a triangulation is achieved by the Delaunay 
triangulation [3]; making use of Fortune's sweepline algorithm [6j, it can be 
computed in 0{n log n). Minimizing the maximum edge can be computed in 
quadratic time, as shown by Edelsbrunner and Tan [H [5]. One of the most 
notorious problems regarding triangulations was finally resolved by Mulzer 
and Rote [9], who proved that finding a triangulations of minimum total 
edge length (a minimum-weight triangulation) is an NP-hard problem. As 
shown by Remy and Steger [11] , there is a PTAS for this problem. The 
maximum-weight triangulation problem has been considered by Qian, and 
Wang [10| . who gave a linear-time approximation scheme for the case of a 
ploint set in convex position, and by Chin, Qian, and Wang [lj, who gave 
a 4.238-approximation algorithm. Computing a MaxMin triangulation has 
been considered by Hu [7], who gave a linear-time algorithm for a convex 
polygon (and thus for a sorted set of points in convex position), and proved 
that the graph version of the problem is NP-hard. Schmidt |12j showed that 
finding a geometric MaxMin triangulation is NP-complete in the presence 
of obstacles, such as inside of a polygon with holes; she also showed that 
computing a MaxMin triangulation for a simple polygon can be solved in 
polynomial time by making use of dynamic programming. 

2 NP-Completeness 
2.1 An Auxiliary Problem 

We start by showing that the following auxiliary problem is NP-complete, 
based on a reduction of Planar 3Sat. 

Problem 1. Covering by Disjoint Segments (CDS) 

Given: A specified set S of line segments ("stabbers") in the Euclidean 
plane, and a subset T of their intersection points ("targets"). 

Wanted: A non-intersecting subset of the stabbers that covers all tar- 
gets. 

This problem is somewhat related to one considered by Megiddo and 
Tamir [8] , who showed that it is NP-hard to compute the minimum number 
of straight lines that are necessary to cover a given set of points in the plane. 
Note, however, that CDS considers line segments of finite length that are 
required to be disjoint. 
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Lemma 2. The problem CDS is NP-complete. 

Proof. We give a reduction of Planar 3SAT, which is the subclass of 
3-satisfiability problems for which the variable-clause incidence graph is 
planar — see Figure [I] To this end, we start with an arbitrary 3SAT in- 
stance I in conjunctive normal form, and use it to construct a CDS instance 
(S, T)j that is solvable if and only if / can be satisfied. 
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Figure 1: A planar straight-line drawing of the variable-clause incidence 
graph Gj of the Planar 3SAT instance / = (xi VX2 V £3) A (x2 V £3 VX4) A 

(xi V X2 V X4). 

In the following, we describe details of the construction; an example of 
the resulting arrangement is shown in Figure [2] 

As a first step, the variable-clause incidence graph Gi is embedded into 
the plane, such that all vertices have integer coordinates bounded by 0(n) 
and the resulting edges are represented by line segments; this can be achieved 
in polynomial time by a variety of graph-drawing algorithms. Let War be 
the set of vertices that represent variables, let V^ a be the set of vertices that 
represent clauses, and let -E c j a be the set of line segments that represent 
edges in Gj. Scaling the resulting graph layout by a factor of 0(n 2 ) results 
in an arrangement in which each pair of vertices are at least a distance of 
0(n 2 ) apart. 

In a second step, we replace each variable vertex v x £ War by an even 
cycle of 0(n) intersecting "variable" line segments surrounding the original 
vertex at distance 0(n). Each cycle consists of 5(v) "even" and 5(v) "odd" 
variable segments; their precise location and parity is chosen in a way that 
ensures that an edge v x v c in Gj intersects an odd segment iff variable x 
occurs in clause c in an unnegated fashion, and an edge v x v c in Gi intersects 
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Figure 2: A CDS instance of the 3SAT instance /. Even variable segments 
are shown in bold, odd variable segments are dotted; clause segments are 
shown in thin solid. Dots indicate the target points. Labels indicate the 
original locations of the variable and clause vertices in the drawing of Gj. 



an even segment iff variable x occurs in clause c in a negated fashion. Let 
SVar be the resulting set of variable segments. Moreover, let S c j a be the set 
of line segments obtained by shifting all segments in E ' j by a distance of 
O(n) towards its clause endpoint, such that intersection with an appropriate 
variable segment is maintained. Let S = SVar U <5 c i a - 

Now let T V ar be the set of intersection points of variable segments, let 
T c j a be the set of intersection points of clause segments, and T = T V ar UT c j a . 
We claim: There is a subset C C S that covers all points in T, if and only 
if there is a satisfying truth assignment for /. 

For the "if" part, consider the truth assignment of a variable x. If x is set 
to be true, choose all the odd variable segments for x; if x is set to be false, 
choose all the even variable segments for x. In either case, all intersection 
points of the variable segments for x are covered. Because every clause c 
must have a satisfying literal, picking the segment that connects the clause 
vertex with the corresponding variable does not intersect one of the selected 
variable segments. 
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For the converse "only if" part, consider a set C C S of non-crossing 
segments that covers all points in T. First it is easy to see by induction 
that if for some x, C contains any even segments, it must contain all even 
segments; otherwise, it must contain all odd segments. This induces a truth 
assignment for all variables. Now it is easy to see that a clause vertex can 
only be covered by a clause segment that does not cross a variable segment, 
which implies a satisfying truth assignment. □ 

2.2 Hardness of MaxMin Triangulations 

Before exploiting the construction of Lemma [2] for the main result, we note 
a helpful lemma; the proof is elementary. 

Lemma 3. Let P be a set of points in the plane, and let Pi,Pj € P. A 
triangulation A contains the edge pipj, iff there is no edge in A that separates 
Pi from pj . 

Now we proceed to the main theorem. 

Theorem 4. It is NP-hard to decide whether a set P of n points in the 
plane has a triangulation with smallest edge of length at least c, for some 
positive number c. 

Proof. Consider the arrangment constructed for the proof of Lemma [2] Let 
Q be the set of all end points of segments in S, and T be the set of target 
points. We perturb all points in Q U T by appropriate powers of 1/n, such 
that the only triples of collinear points correspond to segments in S with 
their covered points. For simplicity, we continue to refer to the resulting 
sets as T and S. As a result, we get a set of points and line segments, such 
that any triangle formed by three points or a segment and a point not on 
the segment has smallest height at least some 5 > 0. 

Furthermore, the segments covering a target point t 6 T subdivide its 
neighborhood into four (for t on a variable cycle) or six (for t at a clause 
vertex) sectors. Replace each point t G T by a pair of points t±,t2 at an 
appropriately small distance e « 5 in opposite sectors. See Figure [3j 

We claim: there is a triangulation with shortest edge of length greater 
than e, iff the corresponding CDS instance can be solved. 

For the "if" part, note that by construction, any segment that covers a 
target point t must separate the corresponding point pair {ti,^}- As all 
points in T are covered, all close pairs are separated, and the claim follows 
from Lemma [3j 
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Figure 3: Replacing the target points of the CDS instance by point pairs at 
distance e. 

Conversely, any pair {ti, t%\ at distance e corresponding to a target point 
t must be separated. By construction (and the perturbation argument), 
a line segment t that connects two points p\ and P2 can only get within 
distance <5 of t if £ covers t in the CDS construction. This induces a solution 
to the CDS instance. □ 

We note an important implication of our construction. 

Corollary 5. Let p(x) be some polynomial. Then the existence of a polyno- 
mial-time algorithm that yields a p(n)- approximation for the problem of find- 
ing a MaxMin length triangulation implies P=NP. 

Proof. In the proof of Theorem [4j choose e small enough that 5/s > p(n). 
Then a p(n)-approximation requires finding a triangulation in which all the 
e-edges are intersected. □ 

3 Conclusions 

Even though our proof implies that finding an approximately maxmin tri- 
angulation in deterministically polynomial time is a hopeless task, there are 
a number of interesting issues that remain. These include practically useful 
methods for constructing exact or approximately optimal solution, as well 
as positive results for special cases and variations. 
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